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Abstract. Dust has been detected in the recurrent nova RS Ophiuchi on 
several occassions. I model the historical mid-infrared photometry and a recent 
Spitzer Space Telescope spectrum taken only half a year after the 2006 eruption. 
The dust envelope is little affected by the eruptions. I show evidence that the 
eruptions and possibly the red giant wind of RS Oph may sculpt the interstellar 
medium, and show similar evidence for the recurrent dwarf nova TPyxidis. 



1. Dusty stellar winds 

1.1. Asymptotic Giant Branch stars and red supergiants 

Stars with an initial mass between ~ 1 and ~ 40 Mq become hydrogen/helium- 
shell-burning Asymptotic Giant Branch (AGB) stars or core-helium-burning red 
supergiants (RSG). These phases are characterised by cool, molecular atmo- 
spheres giving rise to M spectral types (S or C for some chemically peculiar 
AGB stars), strong radial pulsations of this atmosphere on timescales of a year 
or more, and a high luminosity (L > 2,000 L ). T he pulsation may lift th e 
atmosphere high enough such that dust can condense ( Bowen Willson1ll99ll ). 



Radiation pressure on the grains then drives a dust wind which, via collisions 
with molecular hydrogen, drags the gas along with it. At a typical wind speed 
Voo ~ 5 to 30 km s -1 these stars lose mass at rates from M ~ 10 -7 to over 10~ 4 
M yr _1 ( van Loon et al. 19991 ). In AGB stars this leads to the removal of the 



mantle and the premature death of the star, leaving the truncated core behind 
as a cooling white dwarf. The more massive RSGs either explode or evolve back 
to higher surface temperatures, and it is not yet clear how much mass they will 
have shed during the preceding RSG stage. 

Mass-loss rates of cool, luminous stars are most easily estimated from the 
reprocessed radiation emitted by the dust grains at infrared (IR) wavelengths, 
which is particularly conspicuous in the A ~ 10 to 40 /iin region. The main 
problem with this method is that the dust con stitutes only a minor fraction of 
the total mass in the wind, typically ~ 1:200 ( Knapp 19851 ) but poorly known 



for all but the dustiest stars. Carbon monoxide has strong rotational transitions 
at mm wavelengths which can be detected in relatively nearby stars, but despite 
being the most abundant molecule after H2 this too is a trace species of which the 
mass fraction is uncertain — and modelling the CO line requires knowledge of 
the temperature profile throughout the envelope which depends, amongst other 
things, on the dust content. 
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1.2. First ascent Red Giant Branch stars 

Low-mass stars (-Mhutiai < 2 Mq) evolve along a first ascent Red Giant Branch 
(RGB) as hydrogen-shell burning stars with an inert helium core. They reach a 
maximum luminosity of only L t ; p ~ 2,000 L Q , and it becomes problematic for 
them to drive a wind. The situation is worsened by the fact that many RGB stars 
do not pulsate strongly and slowly enough to sufficiently increase the scaleheight, 
and the dust condensation is therefore unlikely to reach completeness. With 
a low, uncertain dust:gas mass ratio and possibly higher dilution of the dust 
envelope as it kinematically decouples from the bulk gas, measured mass-loss 
rates will tend to be too low. The threshold below which this happens is not 
well known, partly because of uncertainties in the opacity of the grains as they 
nucleate and grow: iGail fc Sedlmavd (|l987l ) estimate M > 10~ 6 M yr" 1 but 



Netzer & Elitzuil (|l993h place it at M > 10~ 7 M yr" 1 . lJudge fc Stencell (|l99ll ) 



show that RGB mass loss must be driven by another mechanism, but that this 
appears to be similarly efficient as a dust-driven wind. The warmer stars further 
down the RGB have more prominent chromospheres, which may be accompanied 
by the generation of Alfven or acoustic waves that provide a possible alternative 
mechanism for driving a wind. 

Mass-loss rates estimated from the IR emission fo r the brightest, dusti- 
est RGB stars are typically M ~ 10~ 7 to 10 -6 M ([van Loon et al.l 120061 : 



Origlia et aT]|2007l ). Mass-loss rates from M ~ 10~ 9 to 10" 6 M yr" 1 have 



been determined from the blue-displaced cores of strong optical absorption lines 
and emission in some of these lines, but these are very uncertain too because of 
the sensitivity to the precise excitation and ionization conditions in the wind. 
For the most luminous RGB stars the IR and optical estim ates tend to agree 
within an order of magnitude ( McDonald van Loonll2007l ). 



2. The dusty stellar wind of RS Ophiuchi 



RS Ophiuchi has a spectral type of M2 (jFekel et al.l 2000), and I will later show 



that the luminosity is L ~ 444 L . This places it on an isochrone of a star 
with an initial mass of Mj n itiai ~ 1.1 ± 0.1 M (Fig. 1). Although it is a few 
times fainter than the tip-RGB luminosity and not as cool as it gets, mass loss 
is expected at a level of M > 10 -9 M yr -1 and may be accompanied by dust. 

2.1. Mid-IR observations 

RS Oph is usually classified as an S-type system, which by definition is a system 
devoid of dust as the dusty systems are classified as D-type. The latter usually 
contain a Mira-type large-amplitude variable donor star, which is likely to be an 
AGB star. RS Oph is not a Mira-type variable, and definitely on the RGB, but 
it does have circumstellar dust, and the S-type classification is thus misleading. 
Its IR emission at A = 10 //m was detected and attributed to dust for the first 
time in 1970 bv lGeisel et al.l ( 19701 ). which lead them to remark: "The infrared 



emission from novae closely resembles that from many infrared stars in which 
the rate of mass loss is more constant. This suggests that the formation of 
grains is governed by similar processes in all these systems, and that novae are 
key objects in understanding how the dust is formed." 
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Figure 1 . RS Oph compared to isochrones from iGirardi et al.l ((2000) . 



Dust had thus been detected within three years from the eruption in 1967. 
RS Oph underwent two more eruptions, in 1985 and in 2006, and as far as I 
know the system has been observed at mid-IR wavelengths twice more during 
quiescence and again shortly after the 2006 eruption (Fig. 2). Large differences 
in the mid-IR emission might be expected depending on the temporal proximity 
of an eruption, so it came as a surprise to find the contrary. 

RS Oph was detected i n the InfraRed Astronomical Satellite (IRAS) sur- 
vey m 1983 (jSchaeferl fl986). It is listed in the IRAS Point Source Catalogue 
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(Evans et al. 2007) 



(Geisel, Kleinmann & Low 1970) 1996: ISOPHOT-S (Hodge et al. 2004) 



Figure 2. Times of eruptions and mid-infrared observations. 
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Figure 3. IRAS scans through RS Oph at 12 /jm (left) and 25 /urn (right). 



with F12 = 0. 4 Jy and F25 < 0.32 Jy (the upper limits a t 60 and 100 /jm are 
meaningless); iKenvon. Fernandez-Castro. Sz Stencell (jl988l ) estimate F12 = 0.38 
Jy and F 25 = 0.12 Jy. and lHarrison Gehrd (|1994T ) estimate F 12 = 0.41 ± 0.04 
Jy and F25 = 0.14 ± 0.04 Jy. I have re-measured the IRAS scans, and de- 
rive refined^stmia^esjof^F^ F 25 = 0.09:±0.03 Jy (Fig. 
3). IKenvon. Fernandez-Castro. &; Stencell ( 19861 ) estimated F12 = 0.30 Jy and 
F25 = 0.16 Jy from pointed observations carried out with IRAS in September 
1983 (they estimate that free-free emission may contribute ~ 0.03 and ~ 0.04 
Jy at A = 12 and 25 /im, respectively). The IRAS 1 2-//m flux densities are 
in remarkably good agreement with the iGeisel et ah ( 1970) measurement of 
-F10.4 ~ 0.4 Jy (an N-band magnitude of ttin ~ 4.5 to 5.1 mag) a dozen years 
earlier. Then again, there had not been any eruption meanwhile. 

Soon after the end of the IRAS mission, RS Oph put up a show again. 
Curiously, it took more than a decade until new mid-IR measurements were 
made, in March 19 96 with the PHQT- S instrument onboard the Infrared Space 
Observatory (ISO). lHodge et ahl ( 2004! ) analysed the spectrum shortward of A = 
9 /im and found no evidence for dust emission. However, at those wavelengths 
the dust emission is < 0.1 Jy; the uncertainties in zodiacal light correction and 
instrument calibration and sensitivity do not permit the detection of the dust 
emission in the PHOT-S spectrum around 8-11 /im. 

Another decade passed before new mid-IR measurements were made, this 
time with the Spitzer Space Telescope (Spitzer) triggered by the February 2006 
eruption. By September 2006 the free-free and line emission had weakened 
enough to reveal a continuum and spectral features unambi guously identifi- 
able with predominantly silicate-type dust ( Evans et al.l 120071 ) . Remarkably, 



the Spitzer spectrum agrees frightfully well with all previous photometry (Fig. 
4). 
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Figure 4. Photometry and Spitzer spectrum, overlain with a black body 
representing the star and a Dusty model that reproduces the envelope. 



2.2. Modelling the dust envelope 

I used the radiative transfer code Dusty ( Ivezic &: Elitzur 19971 ) to model the 
spectral energy distribution (SED). The main input parameter of interest is the 
optical depth, ry, which is related to the mass- loss rate via the dust:gas mass 
ratio (assumed to be 1:200), luminosity (L = 444 L , measured from the fit 
to the data assuming a distance of 1.5 kpc) and wind speed (computed by the 
model for a dust-driven wind). The second input parameter of interest is the 
temperature at the inner edge of the dust envelope, To- The model SED also 
depends on the dust species used, and on the temperature of the star (for which 
we adopt T = 3, 600 K consistent with its spectral type). A good fit was ob - 
tained using circumstellar silicates from Qssenkopf. Henning. Mathisl ([l992), 
and To = 600 K. The model overestimates the strength of the 10-/um feature, but 
this also depends on the grain prop erties such as size, shape and composition (a 
( Mathis. Rumpl. &; Nordsieckl 119770 size distribution was assumed). 

I thus estimate a mass-loss rate M = 2.3 x 10 -8 M Q yr _1 at a wind speed 
^oo = 8.1 km s _1 . Although silicates can form already at temperatures above 
1,000 K there is a conspicuous absence of such warm dust (T > 600 K). Because 
the dust shell is optically thin the SED does not depend much on the geometry. 
Assuming a spherically symmetric wind, the dust envelope starts at a distance 
of 39 R* (~ 10 AU) from the red giant (the outer radius is not constrained by 
the available data other than that it is S> 10 2 R*); the orbital radius is ~ 1 AU 
(Fig. 5). Although the wind is unlikely to travel at a constant speed Voo before 
dust even forms, it would take 6 years to reach the dust formation zone if it did; 
the binary would have gone round nearly five times and it is thus highly unlikely 
that the binary motion leaves non-axisymmetric imprints on the dust shell. 
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Figure 5. Sketch of the geometry of the RS Oph system, roughly to scale. 



2.3. How typical is the mass-loss rate of RS Oph? 

It has be en suggested that symbiotic stars lose mass more rapidly than sin- 
gle stars (Kenyan 19881 ). However, when comparing RS Oph with other dusty 
stars I find no evidence for a particularly high mass-loss rate. In Fig. 6 the 
differences are plotted between the measured mass-loss rate and the one pre- 
dicted by the formula for dust-driven winds of oxygen-rich AGB stars and RSGs 
(jvan Loon et all 12005). RSOph sits at the warm end of the diagram, where 
many stars appear to have lower mass-loss rates or a lower dust content (includ- 
ing the famous RSGs Betelgeuse and /iCep). Yet its mass-loss rate is nearly 
as expected for a well-developed, efficiently dust-driven wind. Although at the 
faint end of the distribution, RS Oph is also in line with w hat is measured for 
red giants in globular clusters ([McDonald van Loonll2007l . their Fig. 13). 

On the other hand, stars that do not pulsate in the fundamental mode, 
like Miras do, have lower mass-loss rates or less complete dust formation (Fig. 
6). Compared to those cooler stars, RS Oph seems to have a rather high mass- 
loss rate. The mass-loss rate of RS Oph may even have been under-estimated, 
if the dust:gas ratio is lower than assumed. Physical mechanisms exist that 
might explain enhan ced mass loss in sym biotic systems, e.g., through equatorial 
collimation or tides (|Mikolaiewskal l2000h . 

The remarkable observation is that the system has been losing mass at a 
(near) constant rate, and that the eruptions have done little to affect the dust 
envelope. This can be understood if the eruptions are highly an-isotropic. 
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Figure 6. RSOph compared with other dusty red giants, and with predictions. 



3. Do recurrent novae affect the interstellar medium? 

Recurrent novae are relatively old systems, and hence they have developed large 
peculiar velocities with respect to the ISM. It is therefore possible that recurrent 
novae carve out a path through the ISM as their winds sweep up ISM material 
and/or the nova eruptions leave a trail of successive shells. 

3.1. The case of RS Oph 

A composite of IRAS images at 12, 25 and 60 /im reveals a dark patch of 1- 
2° diameter around RS Oph (Fig. 7) — the direction of the galactic plane is 
towards the lower-left corner. After scaling to the same mean and spread, a 
median combination of all 12-100 /xm images shows an elongated cut into the 
ISM, with RS Oph near its end and its proper motion pointing at a trajectory 
along it. In the same time it took RS Oph to traverse along the structure (~ 1 
Myr), its wind would have blown a bubble of ~ 1° diameter (Fig. 7). This is 
highly suggestive of RS Oph being responsible for the ISM structure. 
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Figure 7. IRAS composites around RSOph (North is up, East is left); the 
right panel is a median combination of scaled 12, 25, 60 and 100 /im images. 



An Ha map (|Finkbeineil 120031 ) shows what looks like a bubble of ionized 
gas, with RSOph near its western edge (Fig. 8). ICool et ail ( 20051 ) noticed this 
structure when studying the nearby symbiotic binary XX Oph; they estimate a 
radial velocity of v = 2.4 ± 11.8 km s _1 for t he gas and exclu de an association 
with XX Oph. RSOph has v = -40 km s _1 (jFekel et al.ll2000l ) and is therefore 
also not associate d with the extended Ha emiss i on. H owever, combined with 
the dust map from ISchleeel. Finkbeiner. David (|l998ft the ISM structure more 
closely resembles a hole in the dust through which background emission is seen, 
rather than a hot bubble (Fig. 8). Composed of a Be star and a late-M giant, 
XX Oph rather than RS Oph might have created the ISM structure. 



3.2. The case of T Pyx 



IRAS images around the recurrent dwarf nova TPyxidis (jGilmozzi Sz Selvelli 



20071 . and references therein) also show a cavity, with a 2" diameter (Fig. 9). A 
10" diamete r hot shell was linked with the 1944 eruption or an earlier one by 
Shara et al. ( 19891 ); the larger cavity I report here could have originated more 
than 10 4 yr ago, possibly built up during successive eruptions. 



4. Concluding remarks 

Dust in RS Oph has been detected in between several eruptions, and at only 
half a year since the 2006 eruption. The silicate dust species and mass-loss rate 
of M ~ 2 to 3 x 10 -8 M yr -1 are not at odds with the expectations for well- 
developed dust-driven winds of single stars. But because RS Oph has a relatively 
warm photosphere, weak pulsation and not very high luminosity, this might in 
fact suggest that the mass-loss rate is enhanced by the effects of the companion 
white dwarf. The eruptions do not seem to have a dramatic effect on the dust 
envelope. An accretion rate onto the white dwarf surface of M > 10~ 7 M Q yr -1 
would probably require an enhanced flow through the Lagrangian point Li. 
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Figure 8. Combination of the COBE dust map (orange) and SHASSA Ha 
survey (blue), showing RSOph and the nearby symbiotic binary XXOph. 



Considering the pattern in the time i ntervals between historic eruptions, 
including the 1907 eruption ( Schaeferl 120041 ) . one might heuristically expect the 



next eruption to be due around 2015 (Fig. 10). 
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Figure 10. Recent outbursts of RSOph and a heuristic prediction for the next. 



